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An XMM-Newton study of the X-ray binary MXB 1659-298 and 
the discovery of narrow X-ray absorption lines 
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Abstract. We report the discovery of narrow X-ray absorption lines from the low-mass X-ray binary MXB 1659- 
298 during an XMM-Newton observation in 2001 February. The 7.1 hr orbital cycle is clearly evident with narrow 
X-ray eclipses preceded by intense dipping activity. A sinusoid-like OM £?-band modulation with a peak-to-peak 
modulation of ~0.5 magnitude and a minimum coincident with the X-ray eclipse is visible. EPIC and RGS spectra 
reveal the presence of narrow resonant absorption features identified with O viii ls-2p, ls-3p and ls-4p, Ne x ls- 



2p, Fe xxv ls-2p, and Fe xxvi ls-2p transitions, together with a broad Fe emission feature at 6.47_ 



keV. The 



equivalent widths of the Fe absorption features show no obvious dependence on orbital phase, even during dipping 
intervals. The equivalent widths of the other features are consistent with having the same values during persistent 
and dipping intervals. Previously, the only X-ray binaries known to exhibit narrow X-ray absorption lines were two 
superluminal jet sources and it had been suggested that these features are related to the jet formation mechanism. 
This now appears unlikely, and instead their presence may be related to the viewing angle of the system. The 
MXB 1659-298 0.6-12 keV continuum is modeled using absorbed cutoff power-law and blackbody components. 
During dips the blackbody is more strongly absorbed than the power-law. The spectral shape of the 3.6% of 
0.5-10 keV emission that remains during eclipses is consistent with that during non-dipping intervals. 
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1. Introduction 

MXB 1659-298 is a transient b urstin g X-ray source dis- 
covered by SAS-3 (Lewin et al. 1976 ). Observations with 
SAS-3 and HEAO-1 detected a clearly modulated per- 
sistent emission which showed erratic intensity variations 
followed by 15 minute duration X-ray eclipses which re- 



peated periodically (Cominsky & Wood 1984). By com- 
bining data spanning ~2 years an orbital period of 7.11 h 
was obtained (Cominsky & Wood 1989). Following the 



MXB 1659-298 wa s obse rved by RXTE a few days 
later (Wachter et al. 2000) confirming that the source 
had resumed its strong and persistent X-ray emission. 
Contemporaneous observations revealed that the optical 
emission exhibits a 0.2 magnitude eclipse feature coinci- 
dent with the X-ray eclipse. The 0.5-30 keV BeppoSAX 
spectrum obtained in 2000 August was modeled by 
a combination of a soft disk-blackbody and a harder 
Comptonized component (Oosterbroek et al. 2001). The 



discovery using EXOSAT of a second bursting, dipping, 



eclipsing, transient (EXO 0748-676; Parmar et al. |1986|) 
the region of sky containing MXB 1659-298 was observed 
by EXOSAT. However, the source was not present. Indeed, 
until recently MXB 1659-298 has not been observed to be 
active. However, on 1999 April 2, after a 21 year quies- 
cent interval, an 18 mCrab (2-9 keV) source at a position 
consistent with MXB 1659-298 was detected in one of the 
BeppoSAX Wide Field Cameras (in 't Zand et al. |1999[ ). 
Observations a day later confirmed that the V"=18.3 opti- 
cal coun terpart to MXB 1659-298 (V2134 Oph; Doxsey et 
al. 1979] ) had brightened significantly and was exhibiting 
high excitation He n emission, confirming that the source 
was undergoing an outburst (Augusteijn et al. 1999| ). 



residuals to this fit suggest the presence of emission fea- 
tures due to Ne-K/Fe-L and Fe-K in the spectrum. 

Only two Low-Mass X-ray binary (LMXRB) systems 
are known that exhibit both dips and eclipses (about 
10 exhibit only dips). During dips the observed spectral 
changes are inconsistent with simple absorption from cold 
material, as might be expected. Detailed modeling of these 
spectral changes provides a powerful means of studying 
the structure and location of the emitting and absorbing 
regions in LMXRB (e.g., Parmar et al 
Balucihska-Church 1995; Church et al. 



1986| ; Church & 
1997T This mod- 
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eling has revealed that at least two emission components 
are required which appear to undergo different amounts 
of absorption during dips. In the case of EXO 0748-676, 
(the other eclipsing and dipping LMXRB), recent XMM- 
Newton Reflection Grating Spectrometer (RGS) obser- 
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Table 1. XMM-Newton instrument configurations 



Instrument 


Energy 


Net 


Mode 




Range 


Exposure 






(keV) 


(ks) 




EPIC PN 


0.2-15 


21.5 


Small window 


EPIC MOS1 


0.2-15 


30.0 


Timing 


EPIC MOS2 


0.2-15 


30.0 


Full window 


RGS 


0.35-2.0 


31.1 


Spectroscopy + Q 


OM 




20 x 1 


B filter 



vations have revealed the presence of discrete structures 
due to absorptio n and emission from ionized Ne, O and 
N (Cottam et al. 2001). Simultaneous European Photon 
Imaging Camera (EPIC) data revealed that the 0.2- 
10 keV spectra can be fit with a two component model con- 
sisting of a central Comptonized component and a more 
extended thermal halo (Bonnet-Bidaud et al. 2001 ). The 
spectral variations observed during dips were mainly ac- 
counted for by variations in the absorbing column affecting 
the Comptonized component only. 



2. Observations and Data Analysis 

The XMM-Newton Observatory (Jansen et al. |200l| ) in- 
cludes three 1500 cm 2 X-ray telescopes each with an EPIC 
at the focus. Behind two of the telescopes there are RGS 
arrays (den Herder et al. 2001). In addition, a coaligned 



optical/UV monitor telescope (OM, Mason et al. 2001) is 
included as part of the payload. Two of the EPIC imaging 



spectrometers use MOS type CCDs (Turner et al . 2001) 
and one uses a PN type CCD (Striider et al. |200l|) . 

The data analysis reported here used event lists 
produced by the standard XMM-Newton SAS (Science 
Analysis Software, V.5.0.0) EMPROC, EPPROC and 
RGSPROC tasks which were further filtered using 
XMMSELECT V.2.35.4. For the PN only pattern (sin- 
gle pixel) events were selected. Known hot, or flickering, 
pixels and electronic noise were rejected using the SAS and 
the low-energy cutoff was set to 0.2 keV. For the MOS, X- 
ray events corresponding to patterns 0-12 were selected. 

Source counts were extracted from circular regions of 
40" radius centered on MXB 1659-298 for the PN instru- 
ment. Background counts were obtained from similar re- 
gions offset from the source position. RGS source and 
background events were extracted by making spatial and 
energy selections on the event files. Wavelengths were then 
assigned to the dispersion coordinates using the latest cal- 
ibration parameters. In order to ensure applicability of the 
X 2 statistic, extracted spectra were rebinned such that at 
least 20 counts per bin were present and such that the 
energy resolution was not oversampled by more than a 
factor 3. In order to account for systematic effects a 2% 
uncertainty was added quadratically to each spectral bin. 
The photo-ele ctric absorption cross sections of Morrison & 
McCammon (1982), ctmm, are used throughout. All spec- 
tral uncertainties and upper-limits are given at 90% confi- 



dence. Quoted count rates are corrected for instrumental 
deadtime, whereas lightcurves are not. 

XMM-Newton observed the region of sky containing 
MXB 1659-298 as a Target of Opportunity (TOO) be- 
tween 2000 March 21 23:55 and March 22 07:16 UTC. 
Unfortunately, the observation was curtailed for technical 
reasons. In addition, much of the data that were obtained 
was not correctly transmitted to the XMM-Newton data 
center at Vilspa. Currently, these data have not been re- 
covered. An examination of the 6.7 ks of PN exposure 
that is available shows that MXB 1659-298 is clearly de- 
tected undergoing dipping, an eclipse and an interval of 
"persistent" emission. However, an examination of the PN 
spectrum shows that the instrument configuration was ab- 
normal and we await an updated calibration before per- 
forming a spectral analysis of this data. Thus, in this paper 
we present only results from a second observation in 2001 
February. 

Due to the problems encountered in the first obser- 
vation, the region on sky containing MXB 1659-298 was 
re-observed by XMM-Newton between 2001 February 20 
08:28 and 16:52 UTC. Table |l| provides details of the in- 
strument modes used. The thin filter was used with the 
EPIC CCDs. The mean 2-10 keV flux obtained from the 
PN corresponds to a luminosi ty of 1.6 10 37 erg s _1 for a 



distance of 15 kpc (Cominsky |198l[) 



3. EPIC Results 

3.1. X-ray lightcurve 

The 0.5-10 keV MXB 1659-298 background subtracted 
light curve obtained from the PN during the 2001 
February observation is shown in Fig. [l] with a binning 
of 32 s. At the start of the observation the source was 
undergoing an interval of rapid dipping activity, prior to 
an eclipse between 2001 February 20 09:10 and 09:25. 
After the first egress the source displayed a long inter- 
val of "persistent" emission punctuated by an X-ray burst 
at 2001 February 20 10:12 hr. The height of the burst is 
strongly reduced by the binning used in Fig. [|. Fitting the 
burst spectrum with a blackbody, after the subtraction of 
the contribution from the persistent emission before the 
burst, gives a temperature of 1.4+°' i keV and a radius of 
15^5 km. Following an interval of ^6 hr, a second eclipse 
is visible. After the second eclipse there is a short, deep, 
dip just prior to the end of the observation. In the 0.5- 
10 keV energy range the first and second eclipse ingresses 
take 15.2 ± 1.1 s and 25.1 ± 1.1 s (at 68% confidence), 
respectively. Similarly the egresses takes 14.4 ± 0.8 s and 
30.4±1.0 s. The dips are deep with the 2-10 keV flux often 
reaching the eclipses level of emission. The upper panel of 
Fig. [l] shows the PN hardness ratio (counts in the energy 
range 2.5-10 keV divided by those between 0.5-2.5 keV) 
with a binning of 128 s. During the dips, which occur be- 
tween orbital phases ~0. 7-1.0, the source becomes harder 
than during the persistent emission. Such behavior is a 
common feature of LMXRB dips. 
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Fig. 1. The XMM-Newton PN 0.5-10 keV light curve of MXB 1659-298 with a binning of 32 s. Time is hours of 2001 
February 20. Two eclipses, periods of intense dipping activity and a burst at approximately 10.1 hrs are visible. The 
upper panel shows the PN hardness ratio, HR, (counts in the 2.5-10 keV energy range divided by those between 
0.5-2.5 keV) with a binning of 128 s. Orbital phase is indicated above 



MXB 1659-298 EPIC PN Persistent Emission 
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Fig. 2. The 0.6-12 keV PN persistent emission. The left panels show the best-fit count spectrum and the residuals. The 
lowest left panel shows the residuals when the normalizations of the narrow absorption and broad emission features 
are set to zero. The structured residuals between 7.0-10.0 keV, as well as the positive residuals around 2.2 keV are 
consistent with being due to instrumental effects. The right panel shows the best-fit deconvolved photon spectrum 
with the individual components indicated. The cutoff power-law dominates ^5 keV 



3.2. Persistent emission spectrum 

The persistent emission spectrum during the 2001 
February observation was investigated by extracting the 
0.2-12 keV PN events obtained between 2001 February 20 
09:30 and 14:00 UTC, which is an interval free of eclipses 
and obvious dipping activity. Times corresponding to the 
X-ray burst around 2001 February 20 10:12 hrs were also 



excluded. This selection results in a 0.6-12 keV count rate 
of 125.0 ±0.1 s _1 with an exposure time of 12.0 ks. In the 
PN small window mode only a 63 x 64 pixel region of the 
central CCD is read out with a time resolution of 5.7 ms in 
order to reduce the effects of pile-up. This is expected to be 
significant at count rates of ^100 s _1 , c ompa rable to that 
of the persistent emission (Struder et al. |200l|) . The effects 
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of pile-up were investigated by extracting events in an an- 
nulus outside the core of the MXB 1659-298 point spread 
function (PSF) and comparing their spectrum with that 
obtained above. As expected, this showed that the overall 
spectrum is affected by pile-up. To minimize this, events 
from the inner 6" radius core of the PSF were excluded 
from the PN spectrum which is used in the subsequent 
analysis. This leaves a background subtracted count rate 
of 78.90 ±0.08 s" 1 . 

Unfortunately, it is not yet possible to analyze MOS 
Timing Mode data. The MOS2 Full Window mode data is 
strongly affected by pile-up and a spectrum of the events 
in the wings of the PSF was extracted. This had a count 
rate of only 6.16 ± 0.02 s _1 . A comparison with the PN 
spectrum showed good agreement between the two, ex- 
cept at energies ^0.6 keV. A fu rther comparison with the 
RGS continuum (see Sect. 4.1) showed good agreement 
between the MOS and RGS spectra and so PN data be- 
low 0.6 keV were excluded. Following discussion with the 
XMM-Newton Science Operations Center personnel, this 
difference is attributed to residual calibration uncertain- 
ties in the PN small window mode. 

We first tried to fit the persistent 0.6-12 keV PN spec- 
trum of MXB 1659-298 with simple models such as an ab- 
sorbed power-law, cutoff power-law (E~ Q exp — (E c /kT)), 
thermal bremsstrahlung, blackbody, and multicolor disk- 
blackbody (Mitsuda et al. |1984| Makishima et al. |1986| ), 
but none of these provided acceptable fits with values of 
reduced \ 2 > 2 and large and structured residuals in all 
cases. Comptonized emission models, such as the COMPST 



and COMPTT models in xspec (see e.g., Titarchuk 1994), 
were the only single-component models providing reduced 
X 2 values of <2, although significant residuals still re- 
mained at high energies. Indeed, the unacceptable fits can 
be mostly attributed to the presence of narrow negative 
residuals around 6.65 and 6.90 keV, suggesting the pres- 
ence of discrete absorption features. Examination of the 
PN background spectra does not reveal the presence of 
features at these energies suggesting that they are not due 
to improper background subtraction. 

Next, we tried to fit the persistent spectrum with dif- 
ferent models, always including the two Gaussian absorp- 
tion lines at energies of ~6.65 keV and at ~6.90 keV. 
Cutoff power-law, COMPST and COMPTT models were 
tested, giving % 2 values of 724.2 for 253 degrees of free- 
dom (dof), 307.3 for 253 dof and 283.8 for 252 dof, respec- 
tively. In all the above fits, the presence of broad positive 
residuals around 6 keV suggests the presence of an addi- 
tional Gaussian line, but this time in emission. This re- 
duced the x 2 to 302.5 for 250 dof, to 282.3 for 250 dof and 
to 264.3 for 249 dof for the cutoff power-law, the COMPST 
and the COMPTT continuum models, respectively. In all 
three models, this additional Gaussian line is significant 
at a confidence level of >99.99 %. Inspection of the resid- 
uals shows that, in the case of the cutoff power-law fit, 
broad structured residuals ^2 keV are still present, sug- 
gesting the need for an additional continuum component. 
Blackbody and multicolor disk-blackbody were added to 



the above continuum and line models and the fits re- 
peated. The inclusion of these low-energy components to 
the cutoff power-law and in the COMPST models results in 
significantly better fits at a confidence level of >99.99 %. 
In contrast, neither low-energy component improved the 
fit with the COMPTT model and so the best-fit with this 
model gives a \ 2 of 264.3 for 249 dof. The best-fit consists 
of absorbed blackbody and cutoff power-law continuum 
components and is shown in Fig. |^. The best fit parameter 
values, given in Table || consist of a blackbody with tem- 
perature, kTbb, of 1.30^Q Qg keV, and a cutoff power-law 
with a photon index, a, of 1.90 ± 0.02 and a high energy 
cutoff, E c , of 15±f keV for a X 2 of 227.1 for 248 dof. This 
model has been used to model the sp ectra of a number 
of dipping LMXRB (e.g., Church et al. 1997 ), allowing an 
easy comparison with results from other sources. For com- 
parison, if the blackbody is replaced by a disk-blackbody 
with a temperature of 2.12+q'q3 keV and a projected inner 



disk radius of 2.20±g;g| km, 



the x is 252.4 for 248 dof. 

2 



The best-fit 7V H of (3.5±g;f) 10 21 atom cm" 2 is a factor 
^2 higher than the galactic column of 1.8 10 21 atom cm -2 
(Dickey & Lockman |1990p , suggesting the presence of ad- 
ditional absorption. 



3.3. Line spectral features 

The two narrow absorption features found in the persis- 
tent emission PN spectrum have energies of 6.64±0.02 keV 
and 6.90^0 01 keV, similar to the energies of the Fexxv 
and Fexxvi ls-2p transitions which occur at 6.70 keV 
and 6.97 keV. These correspond to the He- and H-likc 
ionization states of Fe and are present over a wide range 
of temperatures and values of the ionization parameter, 
£ (= £/n e r 2 , where L is the ionizing luminosity, n c is the 
electron density and r is the distance from the source). The 
measured offset of ~60 eV is consistent with the current 
energy calibration uncertainty of the PN small window 
mode (F. Haberl, private communication). When the in- 
strumental broadening is removed the upper-limits to the 
line widths correspond to velocities of <3600 km s _1 and 
<2300 km s _1 for the Fexxv and Fexxvi transitions, re- 
spectively. We have searched for absorption edges due to 
Fe XXV and Fe xxvi at 8.83 and 9.28 keV which might 
be expected due to the presence of absorption lines from 
these ions. However, there is no evidence for such features, 
but we caution that this is a spectral region where there 
is significant uncertainty in the PN calibration. 

In order to search for any orbital dependence of the ab- 
sorption features, five PN spectra covering orbital phases 
between 0.05 and 0.96 (where center of eclipse is phase 0.0) 
were extracted. These cover intervals of persistent emis- 
sion and dipping activity (phases 0.65-0.90). The spectra 
were processed in the same way as the persistent emis- 
sion spectrum discussed in Sect. |3.2| . The residuals to the 
persistent emission best- fit (excluding the line features), 
are shown in Fig. ||. The depths of the features appear to 
be variable, but there is no obvious dependence of orbital 
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Table 2. Best- fit spectral parameters for the persistent 
MXB 1659-298 emission with a cutoff power-law, with 
photon index a and cutoff energy E c , and a blackbody, 
with temperature fcTbb an d radius i? b b together with 2 
Gaussian features in absorption and one in emission. En ne 
is the line centroid, a its width, EW the equivalent width, 
and June the intensity. Uncertainties are given at 90% con- 
fidence for one interesting parameter, except for upper- 
limits which are 68% confidence. The blackbody radius 
and normalization have been corrected for the missing flux 
when extraction annuli were used and assume a distance 
of 15 kpc 



Component 


Parameter 




Value 




N H (10 22 atom cm -2 ) 




u - JJ -o.oi 


Cutoff 


a 




1.90 ± 0.02 


power- law 


E c (keV) 




15l 29 




Norm (ph. keV -1 cm" 




0.1681°;°°° 


Blackbody 


fcT bb (keV) 








R bb (km) 




4 O+0.07 


Fe emission 


E linc (keV) 




6-471^ 


feature 


a (keV) 




0.60lg;i 2 




/line (10 -4 ph cm -2 s" 








EW (eV) 




1601™ 


Fe xxv 


Tunc (keV) 




6.64 ±0.02 


absorption 


a (keV) 




<0.063 


feature 


/lino (10 -4 ph cm" 2 s" 




1 O+0.3 
J - 4 — 0.8 




EW (eV) 




°°-20 


Fe xxvi 


Eunc (keV) 




6 90+ 02 


absorption 


a (keV) 




<0.043 


feature 


/line (10~ 4 ph cm -2 s" 




_1 4+o- 3 




EW (eV) 








X 2 /dof 




227.1/248 



4.0 

0.0 
-4.0 
-8.0 

0.0 
-4.0 
-8.0 

0.0 
-4.0 
-8.0 

0.0 
-4.0 
-8.0 

0.0 
-4.0 
-8.0 




0.05-0.20 



iff 



0.20-0.35 



0.35-0.50 




6.0 



0.65-0.95 

i i i i 

6.5 7.0 7.5 

Energy (keV) 



Fig. 3. Fit residuals with respect to the average persistent 
emission best-fit model (see Table |^) when the normaliza- 
tions of the two narrow Fe absorption features seen in the 
PN are set to zero for different orbital phase (indicated) 



phase (see Fig. ^), including during the interval of dipping 
activity. This means that similar amounts of the material 
responsible for this obscuration is present in the line of 
sight throughout most, or all, of the orbital cycle. 

A similar analysis is not possible for the broad emis- 
sion feature at 6.47 keV due to its faintness. Assuming 
that the same energy offset (60 eV) should be applied to 
the emission feature as with the absorption features, gives 
an energy of 6.53l~o; 1 4 keV. Thus, the energy is consis- 
tent with Ka emission from a wide range of ionization 
states between neutral and He-like Fe. Given the broad 
line width (a = 0.60t^;^ keV) it is likely that this feature 
originates from fluorescent excitation of Fe with a range 
of ionization states. 



Oosterbroek et al. (2001) report the possible detection 
of a narrow emission line at 1.06 ±0.08 keV with an EW of 
26 ± 16 eV in the BeppoSAX spectrum of MXB 1659-298 
which they attribute to emission from Fe-L or Ne-K ions. 
A similar feature may be present at 94% confidence in the 
EPIC spectrum and is visible in the residuals of Fig. ||. 
Adding such a narrow feature to the best-fit model given 



in Table || gives an energy of 0.99 ± 0.06 keV and an EW 
of 12+^ eV. 



-20.0 

-40.0 

-60.0 
-20.0 

-40.0 

-60.0 



0.0 



Fe XXV 



Fe XXVI 



0.2 



0.4 0.6 
Phase 



0.8 



1.0 



Fig. 4. The EWs of the two Fe absorption features as a 
function of orbital phase. There is no obvious orbital de- 
pendence, even during the dipping activity at phase ^0.8 
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Fig. 5. The PN eclipse spectrum of MXB 1659-298 and the 
best-fit using the persistent emission spectrum consisting 
of a cutoff power-law and a blackbody. The spectral pa- 
rameters were fixed at the values given in Table |, except 
for the normalizations 



In order to model the spectral changes during the 
MXB 1659-298 dips, we used the "complex continuum" 
approach of Church et al. (1997) which has been used 
to model the spectral changes observed from a number 
of dipping LMXRB. Following Church et al. ( |1997[ ), the 
shapes and normalizations of both continuum components 
were held fixed at their persistent emission values. In 
addition to a galactic column, -/V ga i, which was fixed at 
1.8 10 21 atom cm~ 2 (Dickey & Lockman |l99"o| ) , the cutoff 
power-law was obscured by an additional partial covering 
component during the dips. The three line features seen 
in the persistent emission by the PN were also included in 
the model. This gave unacceptable fits with \ 2 values of 
310.3 for 204 dof (shallow dips), 624.7 for 226 dof (medium 
dips) and 564.7 for 191 dof (deep dips). Inspection of the 
residuals shows that a component with a shape similar 
to that of the blackbody is still visible at low-energies, so 
partial covering was added to this component as well. The 
intensity during the dips is then modeled as: 



3.4. Eclipse spectrum 

As Fig. Q indicates, there is significant flux remaining dur- 
ing the eclipses. Counts extracted from both the eclipses 
give a background subtracted PN rate of 5.29 ± 0.09 s _1 
for an exposure of 0.98 ks. With this low count rate it is 
not necessary to exclude the core of the PSF, as for the 
persistent emission. 

The eclipse spectrum can be well fit by the same con- 
tinuum model as for the persistent emission with all the 
parameters fixed at the best-fit values given in Table g, 
except for the normalizations of the two continuum com- 
ponents to give a x 2 OI 95.6 for 125 dof (Fig. ||). This 
model gives a 2-10 keV absorption corrected flux of 
2.3 10~ n erg cm~ 2 s _1 which corresponds to a 2-10 keV 
luminosity of 5.8 10 35 erg s _1 at a distance of 15 kpc. 
This is 3.6% of the persistent emission. The 90% confi- 
dence upper-limits to the EWs of the features at 6.47 keV, 
6.64 keV and 6.90 keV are 380 eV, -121 eV and -42 eV, 
respectively, consistent with the persistent values. 

3.5. Dip spectrum 

In order to study the MXB 1659-298 spectrum during 
dipping intervals, 3 intensity selected spectra were ex- 
tracted from intervals that excluded persistent emission 
and eclipses. Intervals corresponding to PN 0.5-10 keV 
count rates (see Fig of <10 s -1 ("deep dips"), 10- 
40 s _1 ("medium dips") and 40-70 s _1 ("shallow dips") 
were extracted, with a time binning of 32 s. The resulting 
background subtracted PN count rates are 7.21±0.09 s _1 , 
14.7 ± 0.1 s -1 , and 50.3 ± 0.3 s _1 with exposure times of 
1.0 ks, 1.5 ks and 0.5 ks, respectively. The core of the PSF 
was only excluded from the medium and shallow dip spec- 
tra, since the low count rate allowed the use of the entire 
PSF during the deep dips without significant pile-up being 
present. 



e -<TMMN sal [[/ bb e— MM^ bb + (1 _ /bb)] /bb+ 

[/c P ie- ffMMA,H -> + (1 - /cpOKcpl + I Gl + I G2 + I G3 ], 

where Nn bh and -/Vh cp1 are the additional absorptions for 
the continuum spectral components and /bb and / cp i are 
their covering fractions (0 < / < 1); /bb, ^cpl and/Gi are 
the normalizations of the blackbody, cutoff power-law and 
the three Gaussian lines, respectively. Thus, the only free 
parameters in the fitting are A^H bb ) ^H cpJ , /bb and / cp i. 
The fits are significantly better with x 2 values of 276.4 for 
203 dof, 241.6 for 225 dof, and 262.2 for 190 dof for the 
shallow, medium and deep dipping spectra, respectively 
(Table ||). The blackbody suffers significantly more ab- 
sorption than the cutoff power-law (~1.5 10 24 atom cm~ 2 
compared to ^5 10 23 atom cm -2 during the deepest dips), 
but has a smaller covering fraction. Whilst the fit to the 
medium dipping spectrum is reasonably good, the fits to 
the other two dip spectra are formally unacceptable and 
so the results should be treated with caution. 

Inspection of the fit residuals revealed structured pos- 
itive residuals around 6-7 keV. Better fits are obtained 
if a partial covering absorber, with an additional column 
of Nh au and partial covering fraction / gau affects the Fe 
emission. If these parameters are set to the values derived 
for the cutoff power-law, then the x 2 values reduce signifi- 
cantly in the case of the medium dipping interval. Smaller 
reductions in x 2 are evident for the deep and shallow dip 
spectra. This suggests that the iron line does indeed suffer 
extra absorption during dipping intervals, but that we are 
unable to conclude whether the absorbing material has 
similar properties as that responsible for the absorption 
of the blackbody or cutoff power-law components. 

4. RGS Results 
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Table 4. Line features visible in the RGS persistent emission spectrum. EW is the line equivalent width. The wave- 
length uncertainties do not include a systematic uncertainty of O.OlA. The line widths are la upper-limits. The line 
widths expressed in km s _1 and keV are after removal of the instrumental broadening. Theoretical wavelengths are 
from Mewe et al. (|1985|) 



Ion and 
Transition 


Measured 
Wavelength (A) 


Theoretical 
Wavelength (A) 


EW 
(eV) 


Line intensity 
(ph cm -2 s _1 ) 


Line width (a) 
(A) (km s" 1 ) 


(keV) 


O vm ls-2p 
vm ls-3p 
vm ls-4p 
Ne x ls-2p 


18.95+H? 
16.001™? 
15.21 ±0.01 

in ic+0.02 
1Z.10_ Q1 


18.970 
16.010 
15.180 
12.130 


-2.6 ±0.4 
-1.3 ±0.4 
-1.6 ±0.4 
-2.1 ±0.4 


-(3.9 ±0.6) 10~ 5 
-(3.2 ± 1.1) 10~ 5 
-(4.5 ± 1.1) 10" 5 
-(1.0 ±0.2) 10" 4 


<0.045 
<0.070 
<0.058 
<0.049 


<600 
<1200 
<1000 
<1000 


<20 
<80 
<60 
<70 



°< 

T 
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Fig. 7. The first order RGS spectra of MXB 1659-298 with a binning of 0.03 A. The upper curve is the persistent 
emission spectrum and the lower curve that of the dipping intervals. Line identifications are indicated. An offset of 
0.002 photon cmT 2 s _1 A -1 has been added to the persistent emission spectrum for clarity 



Table 3. Best- fit spectral parameters for the 3 dipping 
intervals. The continuum model used is a cutoff power- 
law and a blackbody to which three Gaussian lines are 
added. Only the parameters for the partial covering (PC- 
fabs) models used to represent the additional absorption 
during dips are given. Nh is in units of 10 22 atom cm~ 2 



Parameter 


Shallow 


Dip Spectrum 
Medium 


Deep 


/bb 

^H cpl 

/cpl 

X 2 /dof 


25±^° 
0.74 ±0.09 
1 4 +o.i 

0.75 ±0.01 
276.4/203 


0.67 ±0.02 
14 ±1 

0.9681°:°°; 

241.6/225 


15713 
0.88 ±0.01 
52lt 
0.960 ± 0.002 
262.2/190 



4.1. Absorption lines 

The first order RGS1 and RGS2 summed spectrum corre- 
sponding to the same persistent emission interval as used 



for the EPIC PN was extracted and is displayed in Fig. 0. 
The spectrum is dominated by a smooth continuum with 
a number of low-intensity features evident. There is no 
evidence for the He-like triplet of O vn or O vm l s-2p i n 
emission as seen from EXO 0748-676 (Cottam et al. 2001 ). 
Instead, we detect 4 narrow absorption lines which we 
identify with NeX ls-2p transition at 12.130 A, Ovm ls- 
2p at 18.970 A, Ovm ls-3p at 16.010 A and Ovm ls- 
4p transitions at 15.180 A. These features are shown in 
Figs. ^ and || and their properties summarized in Table |[ 
In addition, there are unidentified features near 23A which 
may be instrumental in origin. The two ls-2p features are 
clearly evident, while the other two O vm transitions are 
less clear. Their measured wavelengths are consistent with 
the theoretical values when the uncertainty in the wave- 
length scale of the RGS (currently 0.008 A, den Herder 
et al. 2001) is included (see Table ^). All the features are 
unresolved and their widths are compatible with the RGS 
line spread function of ^0.06 A. When the instrumen- 
tal broadening has been subtracted, the upper-limits to 
the widths correspond to velocities of <600 km s _1 and 
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Fig. 6. EPIC PN 0.6-12 keV dipping photon spectra de- 
convolved using the partial covering model with the pa- 
rameters given in Table ^. The contributions of the cutoff 
power-law (which dominates at high energies), blackbody 
and line features are shown separately 

<1000 km s^ 1 for the O and Ne ls-2p transitions, respec- 
tively. Fig. H shows the spectral regions around the lines 
in more detail. 

RGS spectra during dipping intervals (orbital phases 
0.65-0.95) were extracted and processed in the same way 
as for the persistent emission spectra. The exposure time 
is 13 ks and the mean count rate 0.72 s _1 . Fig. [?] shows 
the dip spectra together with those obtained during the 
persistent intervals. While the persistent spectra clearly 
show the two lines in absorption, these are not so evident 
in the dip spectra. However, fitting narrow Gaussians to 
the dip spectra gives upper-limits to the EWs of 1.4 ± 0.5 
and 1.5 ± 1.0 eV to any features at the wavelengths of the 
O vni and Ne x ls-2p transitions. Thus, we cannot exclude 
the presence of these features during dipping intervals with 
the same EWs as during the persistent emission. 
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NeX 
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Fig. 8. Persistent emission RGS spectra in the regions 
around the narrow absorption features with a binning of 
0.03 A in units of photon cm -2 s -1 A -1 . Theoretical wave- 
lengths are indicated 



5. OM Results 

A total of 20 B filter OM images were obtained. The coun- 
terpart to MXB 1659-298 (V2134 Oph) is clearly recogniz- 
able in the images and is significantly brighter than in the 
/-band image of Wachter et al. ( 2000[ ). The magnitudes 
of V2134 Oph were obtained by integrating the counts in 
a circular region with a radius of 4" (8 pixels). This is 
slightly smaller than the recommended value of 12 pixels 
(Gondoin 2001), but was chosen because of crowding in 
the field. The background was determined from a region 
0.'45 East of the source with a radius of 16 pixels. The 
background subtracted counts were converted to magni- 
tudes using Star 1 which has a B-magnitude of 18.1 in 
the USNO A2.0 catalog. Due to the relatively poor OM 
spatial resolution and the large difference in brightness, 
the contribution of the faint star close to, and partially 
resolved from MXB 1659-298, is included in the magni- 
tude determinations. The uncertainty in magnitude was 
estimated from a comparison of the derived Star 1 and 2 
magnitudes to be ~2% or 0.02 magnitudes. 

Fig. ^ shows the lightcurve together with the rebinned 
0.5-10 keV PN lightcurve. A sinusoid-like optical modu- 
lation is present, with a narrow minimum coincident with 
the X-ray eclipse superposed. The peak-to-peak amplitude 
is ~0.5 magnitudes and a mean value ^18.5 magnitudes. 
The overall shape of the modu lation is consistent with 
that reported in Wachter et al. ( 2000 ). The broad modu- 
lation of the optical lightcurve with a minimum coincident 
with the X-ray eclipse is consistent with the obscuration 
of an extended object, such as an accretion disk, centered 
on the compact object. 
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Fig. 9. EPIC PN (0.5-10 keV) and OM [B filter) 
lightcurves of MXB 1659-298. Time is hours of 2001 
February 20. An X-ray eclipse is visible at 9:20 hr 



6. Discussion 

We report the results of an XMM-Newton observation of 
the dipping, eclipsing, bursting LMXRB MXB 1659-298 in 
2001 February. EPIC and RGS spectra of the persistent 
emission reveal the presence of narrow resonant absorp- 
tion features identified with O VIII ls-2p, ls-3p and ls-4p, 
Ne X ls-2p, Fe xxv ls-2p, and Fe xxvi ls-2p transitions, 
together with a broad Fe emission feature at 6.47^q keV. 
This is in striking contrast to EXO 0748-676 where RGS 
observations revealed the presence of strong emission fea- 
tures from O VII and O vin as well as weaker emission 
lines from Ne ix, Ne x and N vn and photo-electric 
ab sorpti on edges of both O VII and O VIII (Cottam et 
al. 200l| ). The EWs of the Fe absorption features observed 
from MXB 1659-298 show no obvious dependence on or- 
bital phase, even during dipping intervals. The EWs of 
the O VIII and Ne X ls-2p features are consistent with 
having the same values during persistent and dipping in- 
tervals. The absorption features are narrow with upper- 
limits to any widths as low as <600 km s~ x (O vin, ls- 
2p). The overall shape of the 0.6-12 keV continuum is 
similar to many other LMXRBs and may be modeled us- 
ing blackbody and cutoff power-law components. During 
dips the blackbody suffers significantly more absorption 

^1.5 10 24 atom cm~ 2 com- 



than the cutoff power-law 
pared to 

Similar behavior has been observed from other LMXRBs 



^5 10 23 atom cm 2 during the deepest dips). 



(e.g., XB 1916-053, Church et al. |1997| ), but is in contrast 
to recent XMM-Newto n EP IC results from EXO 0748-676 
(Bonnet-Bidaud et al. |2001 ) and BeppoSAX results from 
Oosterbroek et al. [2001 where the hard components in 
these sources are seen to undergo strong absorption during 
dips. The BeppoSAX result is especially intriguing since 
it is from the same source, with the same quoted 2-10 keV 
luminosity (1.6 10 37 erg s -1 ), and from an observation per- 
formed only 6 months earlier during the same outburst. 
This change in dipping properties, if confirmed, suggests 



21 



20 r 



o 



8 r 
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Log kT (keV) 

Fig. 10. Ion column densities derived from the curve of 
growth analysis for the O vm ls-2p, Ne x, Fe xxv and 
Fe XXVI absorption lines. The curves have been offset 
slightly along the temperature axis for clarity 



that the location, size and structure of the emitting re- 
gions in LMXRBs can change on a timescale of months 
without any apparent change in the X-ray luminosity. 

Narrow X-ray absorption lines have been detected 
from the superluminal jet sources GRO J1655-40 (Ueda 
et al. |1998|; Yamaoka et al. |200ip 
Lee et al 



(Kotani et al. 2000 



and GRS 1915+105 
20011). ASCA observations 



of GRO J 1655-40 revealed the presence of absorption fea- 
tures due to Fe xxv and Fe xxvi which, like those ob- 
served here, do not show any obvious dependence of their 
EWs on orbital phase. GRO J1655-40 has been ob served 
to undergo deep absorption dips (Kuulkers et al. 1998) 



consistent with observing the source at a high inclination 
of 60°-75° (e.g., Frank et al. |1987| ). ASCA observations 
of GRS 1915+105 revealed, in addition, absorption fea- 
tures due to Ca xx, Ni xxvn and Ni xxviii. A recent 
Chandra HETGS observation of this source showed broad 
(~1000 km s _1 ) absorption features from Fe xxv and 
Fe xxvi. Until now, it was possible that these absorp- 
tion features are peculiar to superluminal jet sources and 
related in some way to the jet formation mechanism. With 
the discovery of absorption features from the LMXRB 
MXB 1659-298 this appears not to be the case, and as 
proposed by Kotani et al. ( [2000 ) ionized absorption fea- 
tures may be common characteristics of accreting systems 
which are viewed close to edge-on. The discovery of these 
features from the two superluminal sources may simply 
be a selection effect due to their high fluxes (~1 Crab) 
compared to other systems. 

In the lik ely h ypothesis of photo-ionization (e.g., 
Bautista et al. 1998j ) , the presence of O vm and Ne x ions 
is indicative of conditions where £ is ~100 erg cm s _1 , 
whereas Fe xxv and Fe XXVI are present in re gions where 
£ £ 1000 erg cm s" 1 (Kallman & McCray [l982|) . This 
broad range of £ suggests either that the absorbing ma- 
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terial is present over a wide range of distances from the 
central source, or has a large range of densities. The ab- 
sorption features observed from MXB 1659-298 are present 
during a wide range of orbital phases. This excludes 
the possibility that the absorbing material is associated 
with that responsible for the dipping activity centered on 
phase ~0.8. As pointed out by Ueda et al. (1998) the 
plasma responsible for the absorption features must be 
an-isotropically located, or at least illuminated, around 
the central X-ray source, otherwise emission from outside 
the line of sight would completely fill in the absorption 
features. 

In order to estimate the column densities of O, Ne and 
Fe that produced the observed absorption lines we per- 
formed acurve of growth analysis following that of Kotani 
et al. ( |2000| ) for GRS 1915+105. This allows the EWs of 
the lines to be converted into column densities depending 
on the kinematic temperature of the absorbing material. 
The kinematic temperature includes contributions from 
thermal motions as well as any bulk motions. Absorbing 
columns for the O vin ls-2p, Ne x, Fe xxv and Fe XXVI 
features were derived for a range of assumed temperatures 
(Fig. |Io| ). As expected, a lower kinematic temperature re- 
quires a higher ion column density, and hence a higher hy- 
drogen column density. However the Ah should not exceed 
one Thomson optical depth (1.5 1 24 atom cm~ 2 ) other- 
wise the absorption lines would be strongly diminished. 
We assume abundances of 8.5 10~ 4 , 1.2 1CT 4 and 4.7 1CT 5 
for O, Ne, and Fe, respectively (Anders & Grevesse 1989] ). 
Examination of Fig. O shows that the total Fe column 
density must be 1(R 9 cm -2 for Ah to be less than 
one Thomson depth. This implies that kTp c k5 25 keV. A 
further constraint on the kinematic temperature of the 
absorbing plasma can be obtained from the upper-limits 
to the widths of the lines. In the case of the Fe xxv and 
Fexxvi lines these are not very constraining (<2400 keV 
and <1000 keV, respectively), whilst those fo r the other 
lines, derived using Eq. 1 of Yamaoka et al. fl200l| ), are 
given in Table [J. These show that kTp c is inconsistent 
with the upper-limit obtained from the O vm ls-2p tran- 
sition of <20 keV which implies that the O column density 
is >3 10 17 cm~ 2 . This again suggests that plasma with a 
range of conditions is responsible for the absorption fea- 
tures. Indeed, at a temperature of 25 keV <;50% of the Fe 
will be fully ionized (e.g., Arnaud & Rothenflug 1985 ) as 
will virtually all of the lower Z elements. This probably 
explains why the source is not blanketed by the extremely 
high absorption implied by the Fe lines - most of the mate- 
rial that could photo-electrically absorb at lower energies 
is fully ionized. Indeed, even assuming an extremely high 
kinematic temperature of 1000 keV, the Ah derived from 
the Fe absorption features is 2 10 22 atom cm -2 , a factor of 
10 higher than the excess low-energy neutral absorption 
measured during quiescent in terva ls (Table ||). 

Following Yamaoka et al. (|200l| ) and assuming that the 
plasma is in hydro-dynamical equilibrium in the direction 
vertical to the disk plane, the distance, r, of the absorb- 
ing plasma from the central source can be estimated from 



its thermal temperature using (h/r) 2 x(GMmu/r)=kT i 



th 



(e.g., Shakura & Sunyaev 1972). Here M is the neutron 
star mass, G the gravitational constant, mh the mass of 
a hydrogen atom and h the scale height of the plasma. 
Assuming that h/r = tan(90°— i), where i is the binary 



inclination angle 



80°) gives r Fc < 2.410 s cm for 



kTp c ,> 25 keV. Similarly, the lower limits to r derived 
from the temperatures given in Table | are ro <; 3 10 cm 
and tno i3 9 10 7 cm. For a companion star mass of 0.5 M Q 
the binary separation is 2 10 11 cm and the radius of the 
accretion disk (assumed to be = 0.8 • i?L, where i?L is the 
Roche Lobe radius of the neutron star) is 5 10 10 cm. Thus, 
the material responsible for the Fe absorption must be lo- 
cated close to the central source, at a radius ~200 times 
less than that of the accretion disk, while the material re- 
sponsible for the O and Ne absorption is located at larger 
radii, consistent with the lower value of the ionization pa- 
rameter for these elements. Using the ionization param- 
eter values (~100 for oxygen and neon; >1000 for iron), 
and the distance from the central source, r, calculated be- 
fore, it is possible to roughly estimate the electron density, 
n c : we find n o ~2xl0 18 cm" 3 from the O absorption lines, 
n o ~2xl0 19 cm -3 from Ne and n e ~3xl0 17 cm -3 from Fe. 
Assuming that Nn=n e Ar, where Ar is the thickness of 
the absorbing layer, we derive Ar ^200— 4xl0 4 cm for O, 
Ar ~30-3xl0 4 cm for Ne and Ar <5xl0 6 cm for Fe. 

Since the absorption features occur through a wide 
range of orbital phases, the absorbing material is most 
likely located in a cylindrical geometry with an axis per- 
pendicular to that of the accretion disk. This geometry 
could explain why similar absorption features have not 
been so far observed from other LMXRBs that are viewed 
at angles further from the orbital plane. 
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